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Abstract: For the undamped Kepler potential the lack of precession has historically 
been understood in terms of the Runge-Lenz symmetry. For the damped Kepler problem 
this result may be understood in terms of the generalization of Poisson structure to 
damped systems suggested recently by Tarasov[T]. In this generalized algebraic structure 
the orbit-averaged Runge-Lenz vector remains a constant in the linearly damped Kepler 
problem to leading order in the damping coefficient. Beyond Kepler, we prove that, 
for any potential proportional to a power of the radius, the orbit shape and precession 
angle remain constant to leading order in the linear friction coefficient. 

1. Introduction 

What happens to orbits subject to linear frictional drag? In typical physical settings, 
such as Rydberg atoms or stellar binaries, the effective frictional forces are nonlinear 
and, typically, lead to the circularization of the orbit. Orbital evolution under linear 
friction is special in that, as we show below, the eccentricity and the apsides do not 
change to leading order in the damping. The purpose of this note is to understand this 
elementary result from the underlying dynamical symmetry of the Kepler problem, thus 
demonstrating the utility of a Hamiltonian notion in its non-Hamiltonian generalization. 

In many astrophysical situations, the secular evolution due to friction of orbits 
in a two body system is towards circular orbits. In an orbit in a central field the 
angular momentum scales with the momentum while the energy generally scales with 
the momentum-squared. Friction, assumed to be spatially isotropic and homogeneous 
but time odd, typically scales the momentum. This means generally that the resultant 
secular evolution in central force systems is that in which the energy is minimized at 
fixed angular momentum. This is clearly the circular orbit. The velocity dependence 
of the frictional force is quite relevant, in particular as referenced against the velocity 
dispersion of the (undamped) motion in that central potential. Clearly, under the action 
of such dissipative forces, a consequence of symmetry is that the fiow in orbital shape 
(not size!) has two fixed points, circular orbits and strictly radial (infall) orbits. 

Few physical problems have received more scrutiny than bounded orbits in the two- 
and few- body system. Among these, the two-body Kepler problem is arguably the most 
experimentally relevant and best studied example, having been illuminated by intense 
theoretical inquiry spanning hundreds of years leading to important insights even in 
relatively recent times [21 [3l 11 El El El S 1 M M [HI [IS [Til [151 [161 [HI 
We do not present a systematic review or histiography of this celebrated problem 
(though we thankfully acknowledge also [23l[2ll[25l[28l[26l[271[29l[30] which we have 
found quite useful for our study). We do not aim to contribute to the vast literature 
on astrophysically and microphysically relevant models of friction in orbital problems 
(though the interested reader may find references [3ll[32l[33l[3l[35l[36l[371[38l[39lllQl 
HH HSl Us] a useful launching point for such review) . 

Instead, the our purpose here is to accomodate from the dynamical symmetry group 
point-of-view the result that linear frictional damping (to leading order) preserves the 
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orbit's shape. Although Hamiltonian systems may lose dynamical symmetry completely 
when dissipative forces are included, it can be shown that some structure may remain 
under a modified symplectic form. After a brief introduction to the method by which 
Tarasov extends symplectic structure of Hamiltonian mechanics to dissipative systems, 
we apply it to the determination of the time averages of dynamical quantities. Damping 
invariably introduces new dynamical timescales and the time averaging we implement is 
over times short compared with these timescales (but still long compared with the orbital 
timescales in the undamped problem). Tarasov's construction reveals the relevance of 
the dynamical symmetry algebra to the damped Kepler problem. 

We then compare this aproach to the classic "variations of constants" method of 
orbit parameter evolution by describing an improvement that follows from our study. 
The elementary method can be generalized to non-Kepler homogeneous potentials and 
also determines orbital shape evolution for linearly damped Kepler orbits beyond leading 
order. 

2. Dynamical Symmetry and Tarasov's Construction 

In the undamped Kepler problem the lack of precession is generally understood as a 
consequence of a dynamical symmetry, the celebrated so(4) symmetry formed from the 

— * 

two commuting so(3), one from the angular momentum L = f x p the other from 
the Runge-Lenz vector, S = L x p + ([HllTllH]) being the maximal set of local, 

algebraically independent operators that commute with the Hamiltonian, H = ^-|-\/(r), 
with V{r) = kr"" for /c < and a = —1. (though see [22] for a more precise and general 
statement of the connection between algebra and orbits in a central field). 

{Li,Lj} = 2eijkLk {Li,Sj} = 2eijkSk {Si,Sj} = —2HeijkSk (2.1) 

The length of S is proportional to the eccentricity (and points along the semi-major 
axis of the orbit, in the direction to the periastron from the focus). Defining L and 
S has utility beyond their being constants in the 2-body Kepler problem, for example, 
parameterizing the secular evolution of orbits under various Hamiltonian perturbations 
[29| m]. This so(4) is one of the maximal compact factor groups of the so(4, 2) (the 
conformal group) extended symmetry formed by L, A, H, the generalization of the 
scaling operator TZ = r ■ p and the Virial operator V = \ — |(9rV^(r) ([T71 [20] ) 

The other central potential posessing an easily recognizable dynamical symmetry 
is the multi-dimensional harmonic oscillator {V as given with /c > 0, a = 2). As is 
well known, the isotropic D-dimensional harmonic oscillator's naive 0{D) symmetry 
is part of a larger U{D) dynamical symmetry. For D = 2 harmonic oscillator, note 
that the U (2) symmetry does enlarge further to a so(3, 2) when including TZ, V and 
their generalization (the virial subalgebra Equation f l4.16p through Equation (14.191) of 
each oscillator alone and closes to a sl{2,R) subgroup of the so(3,2)). Note further 
that it is this later algebra that is isomorphic to the dimensionally reduced so(4, 2) of 
the 3-d Kepler problem, by which we mean the reduction of that algebra to generators 
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associated with the orbital plane only. These considerations can also be understood from 
the KS construction [26l [27] of the Kepler problem, in which a four-dimensional isotropic 
harmonic oscillator is the starting point. In that construction the m(4) = sm(4) x u{1) 
is, of itself, not preserved by the KS construction. Instead, it is the u{2, 2) subgroup of 
the four identical, independent oscillator's sp{8, R) symmetry in which the overall u{l) 
can be isolated as the angular momentum constraint of the KS construction [28]. The 
residual symmetry su{2, 2) ~ so(4, 2) is that of the 3-d Kepler problem. The analytical 
connection between the Kepler problem and the isotropic harmonic oscillator has deep 
historical roots, going back to Newton and Hooke (see [15] and references therein). 
Finally, the geometric construction of the undamped Kepler problem as geodesic flow 
on (spatial) a 3-manifolds of constant curvature relates the so(4) dynamical symmetry 
to the isometry group generated by Killing vector fields on the spatial slice [19], EU [23] . 
These various connections between the Kepler problem and the isotropic harmonic 
oscillator do not lead to a simple structural connection between the associated damped 
problems. 

To leading order in the damping, Kepler orbits subject to linear frictional force do 
not change shape or precess as they decay. It would be satisfying to understand this 
elementary result as a consequence of the preservation of the dynamical algebra under 
linear friction. Although this is reminiscent of the damped N-dimensional harmonic 
oscillator, there is no simple way to relate the damped problems. Since the subgroup 
associated with the shape and precession (through the S) is rank one it is suggestive 
that the entire group structure is preserved to leading order in the linear friction. 

A recent paper by Tarasov[T] suggests a straightforward generalization of the 
Poisson structure to systems with dissipative forces. There are many other approaches 
to addressing structural questions of dissipative systems (for one example, see [16| W7\)- 
We find the approach of [1] to be most useful for addressing questions of the dynamical 
symmetries that survive including dissipation. For completeness we now briefly review 
Tarasov's construction, and apply it to dissipation in the central field problem in the 
following section. 

To preserve as much of the algebraic structure as possible, Tarasov constructs a 
one-parameter family of two forms (that define a generalized Poisson structure) that -in 
a sense- interpolate between different dampings. In the zero damping limit it smoothly 
matches onto the canonical symplectic form. Dimensionally, any damping parameter 
introduces a new time scale into the problem, thus this new interpolating two form must 
also be explicitely time-dependent. Tarasov requires this family of two forms to have 
the following useful properties 

(1) Non-degeneracy: The two form u = ujij{t)dx'^ A dx^ is antisymmetric and non- 
degenerate along the entire flow. The x* are the 27V (local) phase space co-ordinates. 
In positive terms, the inverse uj^^ujjk = Sj exists almost globally 

I Since we do not formulate this entirely in the exterior calculus, we must allow for higher codimension 
singularites that may not be resolvable in the dissipative system. 
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(2) Jacobi Identity: the two form is used to define a new Poisson bracket {A, -B}t = 
uj^^diAdjB that forms an associative algebra. Exphcitely it satisfies. 

{A, {B, C}t}t + {B, {C, A}t}t + {C, {A B}t}t = (2.2) 

Here we use the subscript T' to distinguish this bracket from the Poisson bracket of the 
undamped problem. 

(3) Derivation property of time translation: with respect to this new bracket the 
time derivative of the new Poisson bracket satisfies the derivation property (also called 
the Liebnitz rule) 

^{A B]r = B]r + {A (2.3) 

These requirements are remarkable for several reasons. First, property (1) indicates 
that (2) and (3) are possible. The deeper relevance of property (1) is that we can 
regard the two-form as (essentially a) global metric on the phase space. Property (2) 
indicates local mechanical observables in this 'dissipation deformed' algebra form a lie 
algebra. Property (3) is key to the utility of Tarasov's construction for understanding 
constants of motion in dissipative systems. It stipulates that time development in the 
dissipative system, while no longer just { ,H} (or even { ,H}t), must be compatible 
with the structure of the symplectic algebra in the new bracket and thus the (new) 
bracket of time independent quantities in the dissipative system are themselves time 
independent. Thus, just as in the Hamiltonian case, time independent quantities form 
a closed subalgebra. Note that for a Hamiltonian system property (3) is automatic 
since in that case time translation is an inner automorphism of the symplectic algebra. 
In a dissipative system by contrast the Hamiltonian is no longer the operator of time 
translation, but, if Tarasov's construction can be implemented, time translation is still 
an automorphism of the algebra, and as such may be regarded as an outer automorphism. 
Finally, from property (3) it follows after a brief calculation that the two-form uj must be 
time idependent in the full dissipative system, ^ = 0. In terms of symplectic geometry, 
this is metric compatibility of the dissipative flow. 

To proceed with the construction, consider the general flow = x*(x, t). Again, 
these are not assumed to be Hamiltonian flows. Assuming property (1) and using uj^^ 
to form a bracket {A, B}t = u^WiAdjB , property (2) leads to the condition 

Cu'^^dmCo''' + J'^'dmUj''' + U^'^OmU}'^ = 0. (2.4) 

Total time derivatives and derivatives along phase space directions do not commute in 
the flow, 

[^^,d,]A = -djAd,x'. (2.5) 

Using this and the jacobi identity (12.21) . one sees that property (3) implies a condition 
relating the form uj and the flow x*. 



dt 



diXj - djXi where Xj = ^jk{t)x^ (2.6) 
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Given x\ we proceed by solving fl2.4p for an u'^^ that satisfies fl2.6l) . This completes 
Tarasov's construction. 

We breifiy offer a few further remarks helpful to orient the reader. First, in the more 
familiar context of Hamiltonian flows, there = = H} for a local function H on 
the phase space. For this case we can compute in the Darboux frame and learn that the 
usual symplectic form (automatically satisfying (12. 4p ) is a solution also to (12. 6p since 
the RHS in that case is zero. We recognize the RHS of (12. 6p as exactly the obstruction 
to the flow, X* being Hamiltonian. 

Conformal transformation of the two-form, u = Qu, where is a a scalar function, 
can only relate two solutions of (12. 4p and (12.61) IFF the is a constant of the motion 
^ = 0. For in that case (12.60 indicates that 

-Q^^ij = Xjdi^ - Xidj^ (2.7) 

whereas (12. 4p yields 

ujjkdjfl + UkidjQ + uijdkfi = (2.8) 

so, contracting by a-nd comparing with (12. 7p . we learn that Q must be a constant of 
the motion. Thus, each solution is conformally unique. 

We do not know what conditions on lead to the existence of even one non-singular 
simultaneous solution u of (12. 4p and (12. 6p . Tarasov[T] provides an explicit solution for a 
general Hamiltonian system ammended by a general linear frictional force. The general 
question of the existence of uj{t) for a more general is at this point unclear, but 
beyond the scope of this present effort. 



3. Dynamical Symmetry in a Damped System 

Consider damped orbital motion in a central field; 

X = p (3.1) 

P=-drV--(3{p)p (3.2) 
r 

with r = |x| and V{r) the interparticle potential (throughout we take the reduced mass 
to be normalized to 1). The function /3{p) is some general function parameterizing the 
speed dependence of the damping, and this form of the damping function is the most 
general consistent with isotropy and homogeniety of the damping forces. Note that we 
can understand this set as descending from a limit in which the central mass is very 
much larger than the orbital mass though, as in general, damping does inextricably mix 
the center of mass motion and the relative motion. We call linear damping the choice 
of P constant. 

The Equation ( 12.60 takes the form, 

duJa;p{t) 



dt 



d,j:{ujp^x'') - dp{uJxpX^) = dpUJ^.pi{(3{p)p') (3.3) 
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Again, we do not know if solutions to Equation fl3.3p exist and satisfy Jacobi for every 
choice of However, for P{p) = const, there is a simple solution to Equation (13.31) 

that satisfies Jacobi [1], 

uj^j{t) = e^'ooi, (3.4) 

where u is the usual symplectic form of the undamped Kepler problem. Physically this 
corresponds to the uniform shrinkage of phase space volumes under linear damping. 

Clearly, in going from {, } (Poisson bracket) to the new bracket {, }t the relations 
in Equation (12. ip gain a factor of e~^*. The algebra in the new bracket resulting from 
this simple rescaling is still so(4). The utility of this simple change to the algebra of 
Equation (12. ip (which was for the undamped system) is that it is now compatible with 
the evolution under Equation (13. ip and Equation (13. 2p of the damped system. To see 
this in an example, take the first relation in Equation (12. ip and take the (total) time 
derivative of both sides. Then note '^^^^^^^ ^ = —'2/3{2eijkLk) ^ 2eijkLk', i-^- the usual 
Poisson bracket is no longer compatible with time evolution. Duplicating the previous 
line for {Li, Lj}T = 2e~^^eijkLk one learns that this is compatible with the flow Equation 
(13. ip and Equation (13. 2p . Similarly, one may check that all the brackets in Equation (12.11) 
(after replacing {,} with {, }r) are as well. Also note that {Li,H}T = = {Si,H}T, 
though since brackets with H no longer delineate time evolution, these equations do not 
imply that L and S are constants of the motion in the dissipative system (also clear 
from Equation ( ]3.10p below). 

The critique here is familiar to any attempt to reconcile symplectic structure 
and dissipation; fundamentally. Equation (13. ip and Equation (13. 2p still treat x and 
p differently so that time evolution is no longer an element in the dynamical algebra of 
{,}or{,}^. 

To relax the category of 'constants of the motion' sufficiently for dissipative systems, 
consider to what extent dynamical quantities averaged over some number of orbits 
change on a longer time scale, i.e. on a timescale relevant to the dissipation (note 
1/ Pip) is essentially that timescale). Let <> denote time averages over many orbits, O 
a classical observable, and suppose that u; is a solution to Equation (12.60 and the Jacobi 
identity for the system as in Equation (13. ip and Equation (13. 2p . In general, 

<{0,H}t>=< co^nt) {id.O -i-p- mP)dpO) > (3.5) 
=< u;{t) - ^] + u;^P{t)Pip)pd,0 > (3.6) 

Note that sums are implied in the x,p indices of the uj^^{t), the new symplectic form. 
Above we have used isotropy to rewrite the sum in the first term in terms of the 
(normalized) symplectic trace of uj^P{t) which we denote simply as uj(t). To show one 
intermediate step, integrating by parts and using Equation (13.30 we arrive at 

< {O, H}t >= ^A{utO)+ < u^'^u'^f^d^xp - dpXa + x'diuj^^)0 

+ ^^P(3{p)pd,0-uj^> (3.7) 
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Where A(G) refers simply to the overall change of the quantity G over time T. Finally, 
using the Equation (12.51) and the fact that uj^^ satisfies the Jacobi identity we reduce 
the above to 

< {O, H}t >= ^ A(^0)+ < (a;^^'a,,x" - uj^^^'d^alO + u;'^(3{p)pdpO - > (3.8) 

We now specialize to vector fields of the general form Equation (13. ip and Equation (13.21) 
to find, 

^A{ujO) =< {O, H}t + - u;^'{d,iP{p))pO - mpdpO) > (3.9) 

and so making the RHS zero indicates conserved quantities in the non-Hamiltonian 
system. Again, this last result was derived for general P{p), which assumes only that 
the friction is isotropic and homogeneous. In the linear friction case I3{p) = (3 = const. 
For that case, using O = L/uj'^ in the above equation implies that L/uo are constants of 
the motion in this system. Similarly, taking O = S/uj indicates that AS* is proportional 
to {2(3L/uj)x < up > which, again, is zero to first order in /3. This result then applied 
to the case of bounded Kepler orbits with linear damping indicates that the (orbit- 
averaged) Runge-Lenz vector, and thus the dynamical algebra of the Kepler problem, 
is conserved to leading order in the linear friction coefficient. 

In elementary terms, although angular momentum L and S are constants in the 
Hamiltonian system for V{r) ~ ^ they evolve under linear damping of (13. ip . (13. 2p as, 

L = -(3L S = -2(3L X p. (3.10) 

Note that in the weak damping limit, since L is conserved to (9(/5°), the second equation 
time averages to —2 < f3p > x < L >. Thus, again we learn that if the damping were 
strictly linear (/9 constant) then since < p >= 0, the time average of S is 0, again 
indicating that the eccentricity vector would be conserved to leading order. Note also 
that it is straightforward to integrate the L equation explicitely, finding L = Lqc'^* the 
initial condition Lq being identified now a conserved quantity of the dissipative system. 
We use these results in the next section of this paper to ammend the 'textbook' orbital 
secular evolution equations. 

4. The Damped Kepler Problem 

The previous section suggests that (linear-) damped bounded Kepler orbits shrink but 
retain their aspect ratio and do not precess to leading order in the damping. It is well 
known that superlinear damping does lead to circularization whereas sublinear damping 
leads to infall orbits in the Kepler case. So far this begs the questions of whether this 
generalizes to other central field problems, and, if so, then at what order in the linear 
damping coefficient do orbits undergo shape and precessional change. In this section 
we address both questions, first describing a problem that arises using a time-honored 
pertubative method for treating general perturbing forces in the Kepler problem, and 
second, generalize the result of the preceeding section to a broad class of central field 
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potentials. We then establish in precise terms the fate of Kepler orbits under linear 
damping. 

Consider the usual secular orbital evolution method (called "the variations of 
constants") most common in literature on cellestial mechanics, for example, in [18] 
(Chapter 11 Section 5, pg. 323, though see also the treatments of non-linear friction 
in [ini [501 [51] ) . In the "variations of constants' method, orbital response to an applied 
force F = Rx + NL + BL x x, in the orbit's tilt Q, the orbital plane's axis, i, the 
eccentricity e, the angle of the ascending node u the semi-major axis a and the period 
T = 27r/n (in their notation) evolve following [48], 



dfl nar ,,sinM 

-r = N^- 4.1 

dt VI - e2 Sim ^ ^ 

nar ,^ , ^ ^, 

iVcosM (4.2) 



dt 



duj nc?^\ — . „ ^ r . . .df2 , ^ ^, 

— = [-i?cos^ + 5(l + --)sm^]-cos«— (4.3) 

de 



— = na Vl - e2 [/2 sin ^ + fi(cos^ + cos E)] (4.4) 

CJ.6 



da ^ or„ ae . „ „a 



2 



2na2[i?— ^=sin^ + fi — v^r^l (4.5) 



(4.6) 



and where 

dn 3n da 
dt ~ ~2a dt 

with u = Q \ ijj and for the unperturbed Kepler orbit, = 1 + ecos6', P is the latus 
rectum, and E is the anomaly, i.e. r = P{1 — ecosE). The central angle 6 is found 
via the usual definition of angular momentum. When we specialize these Kepler orbit 
evolution equations to the case of isotropic and homogeneous friction we learn that (see 
[48] . Chapter 11, section 7 but using (3{p)p for T in that reference) 

^ = 2paV(p)p (4.7) 

dcij 2 sin^ ^, , , ^ , 

-rr = m 4.8 

dt e 

and 

^ = 2(cose + e)/5(p) (4.9) 
dt 

We can now specialize further to the marginal case, linear friction P{p) = (3 = const. To 
integrate these equations, note r^^ = L = Lqc"^* and, in terms of the force components, 
= 0, and R = (3{p)p cos v and B = j3{p)p sinv where v is the angle between the radius 
vector and the tangent to the orbit. That angle can be written using the parameteric 




Figure 1. The eccentricity in the actual damped Kepler problem (solid curve) 
compared with the eccentricity from (|4.10p and (|4.1ip (dashed curve) versus time. 
Note that for the later the eccentricity can oscillates through zero and can even, as in 
this case, asymptote to a negative value. 



form of r in terms of the constants of the orbit and the angle 9, (sint; = L/rp and 
cosv = esin 9 /Lp) resulting in a self-contained pair of ODE's in e, 9 and t, 

- = -^(l + 6cose)^ (4.10) 

^ = -2/?(cos ^ + e) (4.11) 
dt 

If we integrate these to leading order in /3 only (by, for example, using the first equation 
to eliminate the time derivative to leading order in j3) we do indeed find that the 
eccentricity is an orbit-averaged constant of the motion. But difficulty arises when we 
try to understand these equations beyond leading order in the damping, as a direct 
numerical integration of the equation set reveals (Figure 1). For a broad set of initial 
angles and small initial eccentricities, the e passes through zero and goes negative. 
For comparison, the eccentricity (i.e. the square root of the length of the S vector) 
computed by numerical integration of the original equations of motion for precisely the 
same mechanical parameters and initial conditions is included on that figure. 

Even if one only wanted to assign importance to the asymptotic change in the 
eccentricity, that asymptotic change from integrating the equation pair f l4.10p , (14.111) 
does not scale correctly with the damping coefficient, as may be checked numerically 
(see jlH] for further admonisions against using the "variations of constants" method 
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over long timescales) . Clearly the "variations of constants" method at higher orders 
in the evolution leads to unphysical results at short and long timescales. The fault is 
traceable to the fact that in higher order there are /3- (the damping coefficient) dependent 
terms in the orbit shape whose contributions are ignored substituting for r using the 
undamped Kepler shape of the ellipse. This substitution is however inectricably part of 
the "variation of constants" method. To further clarify this problem with the "variations 
of constants" method, it is not due to some ambiguity in the eccentricity of a non- 
closed orbit, since eccentricity itself, rendered as the length of the S vector, has a local 
definition. Algebraically, with this definition of the eccentricity, note that e^—2L'^U = k"^ 
in the 1/r potential even under arbitrary damping. A more useful algebraically identical 
form is = 4V^r^ — 2TZ^H, from which, since H is negative for any damping function 
on a bounded orbit, we see immediately that is bounded away from zero. 

We now, in two parts, describe an approach emphasising the secular evolution of 
the dynamical symmetry, that addresses this mismatch with the usual "variation of 
constants" method. For simplicity we focus in the main on potentials with fixed scaling 
wieght a, deined through V{r) = kr"". Orbits in any central potential are characterized 
by a fixed orbital plane and a single dimensionless parameter, the ratio d/c oi the 
perihelion distance d to the aphelion distance c. Let L denote the angular momentum 
so that Veff{r) = ^ + V{r) is the effective potential. Then from Veff{c) = U = Veff{d) 
where U is the total energy for a V{r) of a fixed a, we have. 
If c""*"^ — d'^'^^ — d°' 

J = c^-d^ ^ = ^^-^^^ 

that then can be reduced to a 'dispersion relation' between U and L, 

jja+2 



k2L2. = f(d/') (4.13) 
where / in this case is a monotonic function on [0, 1]. Note also that f{x) = f{l/x). 
We call d/c the aspect ratio of the orbit (related to the eccentricity in the a = —1 case). 
Thus in leading order (only) in the damping we think of the RHS as a function of the 
orbital eccentricity only. In applications, the differential form of (14.131) is particularly 
useful, 

(a + 2)^-2a^]^ = ^Wc) (4.14) 

Equations of this sort are often written down when referring to the secular evolution of 
orbital system (see for example [52] and references therein). As before consider further 
only damping forces that are isotropic and homogeneous; they can be written in the form 
from the previous section, F^rag = —l3{p)p- (to simplify notation we henceforth drop 
the vector symbol over the p denoting by p both \p\ and p, unambigious by context). 
In the limit of weak damping we expect L to be approximately constant so that, time 
averaging, we arrive at < 6L >= — < j3{p) > L to leading order in l3{p). Note also that 
SU = —P{p)p'^ to leading order in /3. 

Note that the time derivative of TZ is the sum of a Poisson bracket with H plus a term 
probprtional to P (see Equation (14.161) ). This is, ^ = 2V + (9(/3) which, averaged over 
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bounded orbits, indicates (the virial theorem) that < V >= 0{f3). Thus for V{r) = kr°' 
this imphes that < p"^ >= a < V > +0{(3) so that < U >= ^ < V > +0(/3), which 
to leading order in (3 in Equation fl4.14p indicates 



Thus restricted to hnear damping (/5 constant), but for any a, the averages in fl4.15l) 
factorize trivially and the aspect ratio is unchanged to leading order under linear 
damping. The Equation 04.151) also indicates that this will, in general, not be the 
case for a velocity dependent damping coefficient. Although for potentials with a fixed 
scaling exponent a there is but one dimensionless parameter (See LHS Equation fl4.13p ). 
the introduction of the damping coefficient (3 introduces new length and time scales, 
indicating that the orbital aspect ratio d/c may be a function of (3 and time. The fact 
that (3 is time odd does apparently not preclude its inclusion to linear order in the 
orbital aspect ratio in general. Thus, we repeat, the conclusion that for any monomial 
potentials linear damping preserves the orbital shape is not a consequence of dimensional 
analysis and discrete symmetries. 

As a final check, note that relation Equation fl4.15p and Equation (13.90 are both 
consistent with the attractors of the secular fiow in the orbital shape. For circular orbits 
is a constant of the motion (again to leading order in (3) and thus the LHS of Equation 
fl4.15p is zero, as expected by symmetry. Note that in contrast to (14.150 in Equation 
(13. 9p the change in the eccentricity is proportional to the eccentricity for any (3{p), and 
since the eccentricity vanishes in this limit its orbit-averaged change by (13. 9p does as 
well. Also the strictly radial infall orbit limit is one in which the inner radius, ci — > 0, 
and so the LHS of Equation (I4.15P being non-zero in this limit looks inconclusive. But, 
by the definition of / via Equation (I4.13P we see that in this limit / — * or / oo 
depending on the sign of a. Thus, by (14.121) . L = and remain zero for any f3{p). In 
Tarasov's formulation, since ( 13. 9p is fully vector covariant for isotropic and homogeneous 
(but otherwise arbitrary /5(p)) the change in the S must be along the vector itself for the 
radial infall case. Furthermore, as indicated in the discussion following (13. 9p . the change 
AS* is linear in L (for any (3{p)) which vanishes in the radial infall case. In summary, 
both prescriptions indicate that circular orbits and radial infall must satisfy < S >= 
for any damping function as expected on the grounds by symmetry. 

Furthermore, it is straightforward to go further and perturbatively show using 
(I4.15P and the equations of motion that /3 = const, is the only shape-preserving damping 
function for the Kepler potential (a = —1). The result is clearly common to all monomial 
central potentials only, as it is straightforward to demonstrate a counterexample in a 
more complicated potential. This is due to the fact that there are no additional length 
scales in the potential and is not the case with other potentials, such as the effective 
potential in General Relativity (where the Schwarzschild radius arises as a second length 
scale in the potential). 

Returning to the rather general statement (13.90 . in the Tarasov formulation, the 
explicit time dependence of a candidate constant of motion O gives a second term which 




(4.15) 



Linear frictional forces cause orbits to neither circularize nor precess 



13 



cancels the last two terms. If the operator has a fixed momentum scaling weight (for 
example, L is weight 1 and S is essentially wieght 2), the last two terms will be of 
that same scaling weight only for the case of linear friction, f3{p) = const. Note that 
this argument does not rule out the existence of additional constants of the motion 
in the dissipative system that scale to zero as one goes to the Hamiltonian limit. 
The argument does, however, certify that in the case of linear friction the original 
Hamiltonian symmetries do survive to leading order in that friction. 

Having shown that linear friction preserves the eccentricity to leading order begs the 
question of what happens in higher order in the damping. In the spirit of the discussion 
after f l4.13p where the Virial played a key role, consider the time evolution of that part 
of the dynamical algebra 

n = {n, H} - /3{p)n = 2V - (3{p)n (4.16) 

V = {V, H} - 2i3{p){2V -H) = -^{drV + ^dr{rdrV))'}Z - 2p{p){2V - H) (4.17) 

H = -2f3{p){2V - H) (4.18) 

and for completeness, we have 

{n,V} = H + V -V + ^drV + ^drirdrV) (4.19) 

To orient the reader to the content of these, first note the Hamiltonian limit {i.e. 
f3{p) hmit) for the Kepler case {a = —1), both < V >^ and < TZ/r^ >^ 
as expected. We thus expect both of these time averages to be atleast proportional to 
some positive power of f3. Now, in the abscence of damping V is time even and TZ is time 
odd. Formally, taking j3 to be time odd preserves this discrete symmetry of the above 
evolution equations. Since we expect the < V > and < 7^ > to be analytic functions of 
(3, it must thus be that < V > vanishes quadratically as /5 — » 0. 

An elementary argument now certifies that the < V > must be nonpositive in the 
damped system. Take /3{p) = (3 a. constant. Consider the radial component of the 
velocity, TZ/r. It must average to zero in the /3 — > limit. Since the damped orbit must 
shrink, we thus expect < TZ/r >~< —C(3 > for some positive quantity C (a function of 
the other orbital parameters, etc.). But now take the evolution equation fl4.16p divide 
by r and time average. Clearly, integrating by parts, < TZ/r >= — < TZ/r"^ r >= 
— < TZ^ /r^ > implying that the time average < 2V/r — (3TZ/r > must also be strictly 
negative. But since < TZ/r > must already be negative, the < V > must also be strictly 
negative in the damped system. 

Specializing to Kepler {a = —1), differentiating = jS'p in time and applying the 
equations of motion of the system with friction, we learn that = — 8/5L^V. Using 
the fact that < V > is negative and order (3'^ and integrating both sides, we learn that 
the asymptotic change in the eccentricity to leading order is positive and also of order 
(note the integral itself scales as 1//5). Furthermore, since these are exact evolution 
equations, we have shown that the integration is well behaved throughout. Thus linear 
friction causes Kepler orbits to become more eccentric by a fixed amount that scales 
with the square of the linear damping coefficient. 
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5. Conclusion 

Typically Hamiltonian symmetries lose thier relevance to the geometry of the trajectories 
when damping forces are added to the Hamiltonian system. If the damping is weak, 
homogeneous and isotropic, then for linear damping in monomial potentials, we have 
shown that orbit-averaged shape is stationary. This can understood most easily through 
Tarasov's generalization of conserved quantities from the Hamiltonian context to the 
non-Hamiltonian setting. This approach also quantifies in precise analytic terms the 
fate and subsequent utility of the dynamical symmetry algebra in the associated non- 
Hamiltonian system. 

There are three main frameworks for understanding orbital motion in a perturbed 
central field. The first is directly from the equations of motion; this admits 
straightforward generalization to the non-Hamiltonian case but somewhat obscures the 
structure and fate of the dynamical symmetry group. The second, namely the KS 
construction, embeds the Kepler orbit problem in the higher dimensional set of harmonic 
oscillators with constraints; this illuminates the dynamical symmetry group but does 
not seem to readily admit a generalization to the non-Hamiltonian system. Lastly, 
the geometrical approach, namely that which associates the Kepler Hamilton equations 
to geodesic flow on manifolds of constant curvature, also illuminates the dynamical 
symmetry group while making the generalization to the non-Hamiltonian case somewhat 
unclear. 

In light of these difficulties, we used Tarasov's framework (and applied to the 
damped central field problem here) for extending Poisson symmetries to dissipative 
systems, emphasising its utility in making crisp connections between dynamics, algebra 
and the geometric character of the solutions. Finally, the dynamical algebra remains 
whole in first order in the linear dissipative system, but flow at higher order is not 
trivial. The secular perturbative method "variation of constants" is not adequate to 
explain this, however an elementary method based on the Virial subalgebra explains the 
change in the shape of kepler orbits in higher order in linear damping. 
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